A new measurement of K^ 4 decay and the s-wave 7T7r-scattering length ag 
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A sample of 4 • 10 5 events from the decay K + — > -k + -k~ e + v e (K e 4) has been collected in experiment 
E865 at the Brookhaven AGS. The analysis of these data yields new measurements of the K e 4 
branching ratio ((4.11 ± 0.01 ± 0.11) ■ 10~ 5 ), the s-wave tttt scattering length (ag = 0.228 ± 0.012 ± 
0.003), and the form factors F, G, and H of the hadronic current and their dependence on the 
invariant ittt mass. 

PACS numbers: 13.20.-v, 13.20.Eb, 13.75.Lb 
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More than thirty years ago it was recognized that 
measurements of the properties of K e 4 decay (K ± — > 
7T + 7r _ e ^ e (z7 e )) would provide important information 
about both the weak and strong interactions. This four- 
body semileptonic decay is particularly interesting be- 
cause the two pions are the only hadrons in the final 
state. It allows studies over a broad kinematic range of 
several form factors describing both the vector and axial 
vector hadronic currents, and uniquely of the low en- 
ergy 7T7T interaction in an environment without presence 
of other hadrons. 

While experimental studies of K e i held promise of 
significant physics insight, the small branching ratio of 
about 0.004% has made precise measurements of the de- 
cay parameters difficult For instance, while the 
possibility of extracting the isospin zero, angular mo- 
mentum zero scattering length ag has long been recog- 
nized ||, it was not until 1977, when the Geneva-Saclay 
experiment Q gathered about 30,000 events, that a mea- 
surement was made of this quantity to 20% accuracy. 

On the theoretical side, chiral QCD perturbation the- 
ory (ChPT) j| makes firm predictions for the scattering 
length. The tree level calculation in ChPT using current 
algebra techniques in the soft-pion limit yields ag = 0.156 
(in units of m„) j§. The one-loop (ag = 0.201 ±0.01 §) 
and two loop calculations (ag = 0.217 0) show a sat- 
isfactory convergence. The most recent calculation ||] 
matches the known chiral perturbation theory represen- 
tation of the 7T7T scattering amplitude to two loops |7]] 
with a phenomenological description that relies on the 
Roy equations |9|,|l0|], resulting in the prediction ag = 
0.220 ±0.005. 

The analysis of the Geneva-Saclay experiment [Q com- 
bined with the Roy equations and the inclusion of periph- 



eral irN — > ttttN data led to the presently accepted value 



of 



= 0.26 ± 0.05 pi. It has been argued, that, if 



the central experimental value ag = 0.26 would be con- 
firmed with a smaller error, such a large value can only 
be explained by a significant reduction of the quark con- 
densate < 0|mu|0 >, as is possible in generalized chiral 
perturbation theory (GChPT) p2| . The quark conden- 
sate is an order parameter measuring the spontaneous 
breaking of chiral symmetry. Vice versa measuring ag 
with higher precision allows to reduce the bounds on this 
parameter 13|. 
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FIG. 1. Plan view of the E865 detector. A K e 4 event is 
superimposed. 

The analysis outlined here is based on data recorded at 
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the Brookhaven AGS, employing the E865 detector. The 
apparatus, described in detail in pi] , is shown in Fig- 
ure Q. The detector resided in a 6 GeV/c unseparated K + 
beam directly downstream of a 5 m long evacuated decay 
volume. A first dipole magnet separated the K + decay 
products by charge. A second dipole magnet sandwiched 
between four proportional wire chambers (P1-P4) served 
as spectrometer. Two gas Cerenkov counters CI and C2, 
filled with CH4 at atmospheric pressure, and an electro- 
magnetic calorimeter distinguished n , and /i^ from e^. 
7T^~ are separated from /i^ by a set of 12 muon chambers. 
Four hodoscopes were added to the detector for trigger 
purposes. In our analysis we determined the K + momen- 
tum using the beam line as a spectrometer, the position 
of the decay vertex and the information from the pixel 
counter installed just upstream of the decay volume. 

The first level trigger selected three charged particle 
tracks based on coincidences between the A-, and D- 
hodoscope and the calorimeter. The second level trig- 
ger indicated the presence of an e + not accompanied by 
an e~ . It required signals in both right side and only 
minimal light in both left side Cerenkov counters. In 
this we discriminated against two of the most common 
background channels: (1) K + — > 7r + 7r + 7r~ (K T ) and (2) 
K + — > 7T + 7r° followed by 7r° — > e + e~7 (Kdai)- 

The offline analysis selected events containing three 
charged tracks with a vertex within the decay volume 
of acceptable quality, a summed momentum of less than 
5.87 GeV/c, and a timing spread between the tracks con- 
sistent with the resolution of 0.5 ns. Even after particle 
identification criteria were applied, the remaining sample 
still contained background events mainly from A" T -decay 
with a misidentification of a n + as an e + and acciden- 
tals. Requiring that the K + reconstructed from the three 
charged daughter particles does not track back to the 
target reduced the background from K T to the level of 
1.3±0.3%, since for K e 4 the undetected neutrino made 
the reconstruction incomplete. The dominating acciden- 
tal background was combination of a tt + tt~ pair from a 
K T decay with an e + from either the beam or a coin- 
cident decay with an e + in its final state. A likelihood 
method was employed to reduce this background to a 
level of 2.4±1.2%. Due to the excellent particle identifi- 
cation capabilities of our detector all other backgrounds 
were negligible. 

After the event selection 406,103 events remained, of 
which we estimate 388270 ± 5025 to be K e 4 events. 

To determine the branching ratio, the form factors 
and other related quantities a Monte Carlo simulation 
is needed. Our code, based on GEANT, takes into ac- 
count the detector geometry as well as the independently 
measured efficiencies of all detector elements. K e 4 de- 
cays are modeled by ChPT on the one loop level jl5Ul6||. 
Radiative corrections are included following Diamant- 
Berger |l7). With this apparatus, we generated 81.6 • 10 6 
K e 4 events, resulting in 2.9 • 10 6 accepted events. The 



agreement between data and Monte Carlo in all control 
variable distributions is very good, as e.g. evidenced by 
the plots shown in Fig. [| 

The K e 4 branching ratio is measured with respect 
to K T decay. K T events were collected in a minimum 
bias prescaled trigger together with K e 4 events. With 
Br(r) = (5.59 ± 0.05)% @, the K ei branching ratio is 
calculated to be 

BR{K ei ) = (4109 ± 8 (stat.) ± 110 (syst.)) • 10" 8 . 

This result agrees well with the average of previous ex- 
periments |§: (3.91 ± 0.17) • 10~ 5 . The systematic 
uncertainties are dominated by the uncertainties in the 
Cerenkov counter efficiencies and background contribu- 
tions. 

The kinematics of K e 4 decay can be fully described by 
five variables jij): (1) s w = M%„, and (2) s e = M% v , the 
invariant mass squared of the dipion and the dilepton, re- 
spectively; (3) 0^ and (4) 8 e , the polar angles of tt + and 
e + in the dipion and dilepton rest frames measured with 
respect to the flight direction of dipion and dilepton in 
the K + rest frame, respectively; (5) 4>, the azimuthal an- 
gle between the dipion and dilepton planes. The FWHM 
resolution of the apparatus for these five variables is esti- 
mated to be: 0.00133 GeV 2 (s^), 0.00361 GeV 2 (s e ), 147 
mrad (6*^), 111 mrad (8 e ), and 404 mrad {(f)). 

The matrix element in terms of the hadronic vector and 
axial vector current contributions and is given by 





280-294 (285.2) MeV 


294-305 (299.5) MeV 


F 


5.832 ± 0.013 ± 0.080 


5.875 ±0.014 ±0.083 


G 


4.703 ± 0.089 ± 0.069 


4.694 ± 0.062 ± 0.067 


H 


-3.74 ±0.80 ±0.18 


-3.50 ±0.52 ±0.19 


8 


-0.016 ± 0.040 ± 0.002 


0.068 ±0.025 ±0.001 


x7NdF 


1.071 


1.080 


Mtttt (A/™) 


305-317 (311.2) MeV 


317-331 (324.0) MeV 


F 


5.963 ± 0.014 ± 0.090 


6.022 ±0.016 ±0.094 


G 


4.772 ± 0.054 ± 0.070 


5.000 ±0.051 ±0.082 


H 


-3.55 ±0.44 ±0.20 


-3.63 ±0.41 ±0.023 


8 


0.134 ±0.019 ±0.002 


0.160 ±0.017 ±0.002 


x7NdF 


1.066 


1.103 




331-350 (340.4) MeV 


> 350 (381.4) MeV 


F 


6.145 ±0.017 ±0.096 


6.196 ±0.020 ±0.083 


G 


5.003 ± 0.049 ± 0.083 


5.105 ±0.050 ±0.074 


H 


-1.70 ±0.41 ±0.024 


-2.23 ± 0.48 ± 0.033 


8 


0.212 ±0.015 ±0.003 


0.284 ±0.014 ±0.003 


x7NdF 


1.093 


1.034 



TABLE I. Form factors and phase shifts 8 = 8q — 6{ for the 
six bins in M W7r . The number of degrees of freedom for each 
fit is 4796. The first uncertainty is statistical, the second 
systematical with dominant contributions from background 
and Cerenkov efficiency. 
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m = ^v:mpMi - isMPeW^ - An , (i) 

A» = FP^ + GQ' 1 + RL^ , V 1 = H^ vpa L v P p Q a , (2) 

where P — p\ ± p 2 , Q = Pi — P2, and L = p e ± p u , and 
Pi, P2, p e , and p v are the four- moment a of the 7r + , 7r~, 
e + , and v e in units of Mk, respectively. The form factors 
F, G, R and H are dimensionless complex functions of 
s n , s e and 9^. The expressions for the decay rate derived 
from this matrix element have been given in ref. pOj . 
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FIG. 2. Phase shift difference 8 as a function of dipion 
mass. The dashed line represents the fit to Eq. ^ for the 
Geneva-Saclay data [0 and the solid line for our data as a 
function of the scattering length a%. 

Amoros and Bijnens recently developed a parameter- 
isation of these form factors,based on a partial wave 
expansion in the variable 9 V [pi| : 



idl 



F=(fs + f' s q 2 + f' s 'q i + feS e ) 

+f P {Q 2 /s,) 1/2 (P-L) cobO* e 
G={g P + g' p q 2 + 9e s e ) e iS * ,H=(h p + ti p a 2 ) e i6 ' , (3) 

where q = {s^ / {AM 2 ) — l) 1 / 2 is the pion momentum in 
7T7T rest frame. The form factor R enters the decay distri- 
bution multiplied by m 2 and can therefore be neglected. 
This parameterisation yields ten new form factors f s , f' s , 
f", fe, fp, 9 P , 9 P , g e , h p , and h' p , which do not depend 
on any kinematic variables, plus the phases <5§ and 5{, 
which are functions of s n . 

The phase shifts can be related to the scattering 
lengths. A recent analysis used the parameterisa- 
tion proposed by Schenk p3]: 



tan 5j = \ 1 



s L — ' 

fc=0 



/4M 2 



V s-si 



(4) 



The Roy equations || are then solved numerically, ex- 
pressing the parameters A\ k and s\ as functions of the 
scattering lengths a% and a 2 ,. The possible values of the 



scattering lengths are restricted to a band in the a" ver- 
sus Oq plane. The centroid of this band, the universal 



curve p3| relates a[j and a, 



at = -0.0849 



0.232 ag- 0.0865 (ag) 2 



(5) 



For the fits we divided our data into six bins in s ff , five 
in s e , ten in cos^, six in cos# e and 16 in (j>. In the x 2 
minimisation procedure, the number of measured events 
in each bin j is compared to the number of expected 
events given by: 



Vi = 



N 



- i> 1- . ' sr MF,G,H) n ™ 



(6) 



where the sum runs over all Monte Carlo events in bin 
j. N K is the number of K + decays derived from the 
number of K T events. N MC is the number of generated 
events. J 5 (F,G,H) MC (= I Q) is the five-dimensional 
phase space density generated at the momentum q = 
q MC with the form factors F, G, and H used to simulate 
the event. J 5 (F,G,H) new is calculated at q MC with F, 
G, H evaluated from the parameters of the fit. Thus, we 
apply the parameters on an event by event basis, and, 
at the same time, we divide out a possible bias caused 
by the matrix element, making the fit independent of the 
ChPT ansatz used to generate the MC. 

For the fit, we have assumed that F, G, and H do not 
depend on s e and that F contributes to s-waves only, 
i.e. f e = g e = f P = 0. Our first set of fits is done 
independently for each bin in s % . The above assumptions 
then leave one parameter each to describe F, G, and H 
aside from the phase difference 5 = 5q — 8±. The results 
are listed in Table |. The centroids of the bin (M^) 
are determined following Lafferty and Wyatt (24|]. If the 
six phase shifts in Table [D are fit using Eq. I and Eq. |, 
one obtains ag = 0.229 ± 0.015 (x 2 /NdF = 4.8/5). The 
resulting curve is shown in Fig. ^|. 



fs 
fs 

g'p 



5.75 ± 0.02 ± 0.08 f' a 1.06 ±0.10 ±0.40 

-0.59 ±0.12 ±0.40 g p 4.66 ± 0.05 ± 0.07 

0.67 ±0.10 ±0.04 h p -2.95 ±0.19 ±0.20 
0.228 ±0.012 ±0.004 



a° 



TABLE II. Form factors and scattering length oq in the 
parameterisation of Eq. | (x 2 /NdF = 30963/28793.) 



fv 
h 

9e 



-4.3 ± 1.3 ± 3.4 
-4.1 ± 1.3 ± 3.1 
0.5 ±4.4 ±11.3 



X = 30952 
X 2 = 30954 

x 2 



30963 



TABLE III. Fit of form factors f p , f e , and g e . The sys- 
tematic uncertainties are dominated by the resolution of the 
neutrino mass squared. (NdF=28792) 
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We have also made a single fit to the entire data sam- 
ple. In this second fit we substituted 5 in Eq. || by the 
expression of Eq. || With the relation between a[J and 
a\ given by Eq. | only f s , f' a , f'J, g p , g' p , hp, and ag then 
remain as free parameters. The results, listed in Table |n| 
are in an excellent agreement with the ones derived in 
the previous paragraph. 

To check the assumption f e — g e = f p = we also 
allowed these form factors to vary, one at the time, in 
our second fit. Table III shows that all three form fac- 
tors are consistent with zero. The quality of the fits is 
demonstrated in Fig. |3|, where the invariant mass (s^) 
and azimuth (</>) distribution from data are compared to 
the reweighted Monte Carlo distributions (Eq. ||) . 
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FIG. 3. Invariant mass and angle distributions describing 
the K e 4 decay. The histograms are the Monte Carlo distribu- 
tions while the markers with the error bars represent the data. 
The dashed histogram indicates the non-A" e 4 background. 

The asymmetry of the <fi distribution is the only ob- 
servable, which directly depends on the phase shifts Q . 
The amplitude of this asymmetry amounts to only 10% of 
the <j) independent part. This explains why the statistical 
error of a® 1S still limited to 5.3%. 

To summarize, experiment E865 has collected a K e A 
event sample more than ten times larger than all previ- 
ous experiments combined. From the model independent 
analysis of this data the momentum dependence of the 
form factors of the hadronic currents as well as 7T7t scat- 
tering phase shifts have been extracted. The form factors 
and phase shifts serve as an important input in the pro- 
gram to determine the couplings of the effective Hamil- 
tonian of chiral QCD perturbation theory at low ener- 
gies [ p5| . From a preliminary communication of these 
results already tight bounds on the value of the quark 
condensate have been extracted |u| . Using the relations 
between a° and a§ given by the Roy equations ]l0[ ], we 
have extracted the most precise value of the itit scatter- 
ing length ag. This value agrees well with predictions 
obtained in the framework of ChPT |8| . 
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